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Spin liquids on a honeycomb lattice: Projective Symmetry Group study of Schwinger 

fermion mean-field theory 
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^Department of Physics, Boston College, Chestnut Hill, MA 02467 
(Dated: July 20, 2011) 

Spin liquids are novel states of matter with fractionalized excitations. A recent numerical study 
of Hubbard model on a honeycomb lattice^ indicates that a gapped spin liquid phase exists close 
to the Mott transition. Using Projective Symmetry Group, we classify all the possible spin liquid 
states by Schwinger fermion mean-field approach. We find there is only one fully gapped spin liquid 
candidate state: "Sublattice Pairing State" that can be realized up to the 3rd neighbor mean-field 
amplitudes, and is in the neighborhood of the Mott transition. We propose this state as the spin 
liquid phase discovered in the numerical work. To understand whether SPS can be realized in the 
Hubbard model, we study the mean-field phase diagram in the Ji — J2 spin-1/2 model and find an 
s-wave pairing state. We argue that s-wave pairing state is not a stable phase and the true ground 
state may be SPS. A scenario of a continuous phase transition from SPS to the semimetal phase is 
proposed. This work also provides guideline for future variational studies of Gutzwiller projected 
wavefunctions. 



I. INTRODUCTION 

Traditional Landau's theorj*^!^ points out that states 
of matter can be classified by their symmetry. And 
the low energy excitations can be understood by either 
bosonic modes or fcrmionic quasiparticlcs, which carry 
integer multiples of the quantum numbers of the fun- 
damental degrees of freedom. Fractional quantum Hall 
liquids (FQHLs) provide a striking counterexample of the 
Landau's paradigm: different FQHLs all have the same 
symmetry, yet they are very different since a quantum 
phase transition is required to go from one liquid to an- 
other. To understand their differences, one has to go 
beyond Landau's paradigm and the concept of topolog- 
ical order was introduced^. The quasiparticle excita- 
tions in FQHLs carry only a fraction of the fundamental 
electric charge. Meanwhile these fractionalized quasipar- 
ticlcs obey neither bosonic nor fermionic statistics and 
are dubbed anyons consequently. 

Can strong interactions lead to similar novel states of 
matter in the absence of magnetic field? After the orig- 
inal proposal of Anderson^, intensive theoretical studies 
have revealed that spin systems can realize such novel 
phases of matter: spin liquids(SL). And a few experimen- 
tal systems have been identified to be likely in spin liquid 
phases^Ti^i. A spin liquid is often defined to be a quan- 
tum phase of spin-1/2 per unit cell that does not break 
translation symmetry. These liquid phases of spins arc 
also distinct from one another by their topological order. 
Although a rigorous theorem is lacking because we are 
still unable to classify all possible topological order, it 
is generally believed that the excitation of a topological 
ordered phase is fractionalizedi^. 

Although theoretical studies have shown that 
spin liquid ground states exist for artificial model 
Hamiltoniansii"— , it remains unclear whether a simple 
or experimentally realizable Hamiltonian can host such 
novel states. Recently a remarkable quantum Monte 



Carlo simulation of Hubbard model on a honeycomb 
lattice^ indicates that a gapped spin disordered ground 
state exists in the neighborhood of the Mott transition. 
Although a honeycomb lattice has two spin-1/2 per unit 
cell, it is impossible to have a band insulator phase 
without breaking lattice symmetry. Therefore this spin 
disordered phase should be topologically ordered and 
have fractionalized excitations. Wc will still term it a 
spin liquid. 

What is the nature of this spin liquid phase? In this pa- 
per we try to propose the candidate states by Schwinger- 
fermion (or slave-boson) mean-field approac h^^d^" — , fol- 
lowing the techniques developed in Ref. [23. Our re- 
sults can be summarized as follows. We search for the 
fully gapped spin liquids which lead us to focus on the 
Z2 mean-field states. We first use Projective Symmetry 
Group (PSG)^ to classify all 128 possible Z2 mean-field 
states that preserve the full lattice symmetry as well as 
time-reversal symmetry. Notice the spin liquid phase in 
the numerical work seems to be connected to the semi- 
metal phase by a second-order phase transition, which 
suggests this state to be in the neighborhood of a uniform 
Resonating- Valence-Bond (u-RVB) state. So we classify 
all the 24 possible Z2 states around the u-RVB states. 
Among these 24 states, we find only 4 states can have 
a full energy gap in the spinon spectrum, while other 
20 states have symmetry protected gapless spinon ex- 
citations. We find that up to 3rd neighbor mean-field 
amplitudes, only one of the four fully gapped Z2 state 
can be realized, and we term it as Sublattice Pairing 
State (SPS). We propose this state to be the spin liq- 
uid state discovered in the numerical study. We also 
study the mean-field phase diagram of the Ji — J2 an- 
tiferromagnetic spin-1/2 model on a honeycomb lattice 
to understand whether SPS can be more favorable than 
the u-RVB state while both states are in the neighbor- 
hood of the Mott transition. We find when J2 > 0.85 Ji 
a spinon gap opens up by s-wave pairing on top of the 
u-RVB state. This s-wave pairing state is not a stable 



phase and is an artifact of the mean-field study where 
gauge dynamics are ignored. On the other hand, the 
proposed SPS Z2 state is continuously connected to the 
s-wave pairing state by making the pairing phase sublat- 
tice dependent. This suggests the ultimate fate of s-wave 
pairing state may be SPS. We propose that a more careful 
projected wavefunction study, which includes the gauge 
fluctuations, may be able to find SPS Z2 state as the 
ground state in the J\ — Ji model. The possible contin- 
uous phase transitions from SPS into semi-metal phase 
are discussed. 
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II. SCHWINGER-FERMION APPROACH AND 

PSG 

In Schwinger-fermion approach, a spin- 1/2 operator at 
site i is represented by: 



Si = ^ha^aHfiP- 



(1) 



A Heisenberg spin Hamiltonian H ~ X]<i7"> Jij^i ' ^j 
is represented as iJ = E<«j> -^•^u(/iL/i"/]/3/j73 + 
^fiafiafjafjis)- Because this representation enlarges the 
Hilbcrt space, states need to be constrained in the phys- 
ical Hilbert space, i.e., one /-fermion per site: 



ft f. 



fiafifS^ap — 0. 



(2) 



Introducing mean-field parameters rjijCafi = —2{fiafjp), 

XijSap = 2(//a/j/3), where eai3 is fully antisymmetric ten- 
sor, after Hubbard-Stratonovich transformation, the La- 
grangian of the spin system can be written as2^ 



i <ij> 

- {iplU.jxl^j + h.c.)] + Yl aUOV'KV* 



(3) 



where two-component fermion notation f/ji = (/i.-j-, /, i )^ 
is introduced for reasons that will be explained shortly. 
Uij is a matrix of mean-field amplitudes: 



[/,;, 




(4) 



00(1) are the local Lagrangian multipliers that enforces 
the constraints Eq.(I5|). 

In terms of ij}, Schwinger-fermion representation has 
an explicit SU{2) gauge redundancy: a transformation 
V-j -^ W^^pi, Uij -^ W,U,jW], W, e SU{2) leaves the ac- 
tion invariant. This redundancy is originated from repre- 
sentation Eq.([T|): this local SU{2) transformation leaves 
the spin operators invarianli^S and thus does not change 
physical Hilbert space. 

One can try to solve Eg.® by mean- field (or saddle- 
point) approximation. At mean-field level, Uij and Og 



FIG. 1: (color online) (a) Honeycomb lattice and generators 
of symmetry group, (b) SPS mean-field ansatz in terms of 
/-fermion. t, A, 9 are real and 6 ^ 0,tt, ±7r/2. The pairing 
phases for A sublattice (blue sofid fine) and B sublattice (red 
dashed line) are opposite. 



are treated as complex numbers, and Oq must be chosen 
such that constraints Eq.® are satisfied at the mean 



field level: {ipli 
written as: 



V. 



0. The mean-field ansatz can be 



Hmf = - ^ iPi'Uiji'j + ^ ipjaQT'^ipi 



(5) 



<u> 



where 



- T U- ■ 



A local SU{2) gauge transforma- 
tion modify Uij —J- WiUijW^ but docs not change the 
physical spin state described by the mean-field ansatz. 
By construction the mean-field amplitudes do not break 
spin rotation symmetry, and the mean field solutions de- 
scribe spin liquid states if translational symmetry is pre- 
served. Different {uij} ansatz may be in different spin 
liquid phases. The mathematical language to classify dif- 
ferent spin liquid phases is PSG—. 

PSG is the manifestation of topological order in the 
Schwinger-fermion representation: spin liquid states de- 
scribed by different PSG's are different phases. It is de- 
fined as the collection of all combinations of symmetry 
group and SU(2) gauge transformations that leave {uij} 
invariant (as Oq are determined self-consistently by {uij}, 
these transformations also leave Cg invariant). The in- 
variancc of a mean-field ansatz {uij} under an clement 
of PSG Girll can be written as 

GuUiiuij}) = {uij}, 

U{{uij}) = {uij = M[/-l(,),,7-l(j)}, 

Gu{{uij}) = {uij = Gu{i)u.ijGu{j)''}, 

Gu{i) e SU{2). (6) 

Here U G SG is an element of symmetry group (SG) of 
the spin liquid state. SG on a honeycomb lattice is gen- 
erated by time reversal T, reflection cr, 7r/3 rotation Cg 
and translations Ti,T2 as illustrated in FIG. [1] (see also 
appendix [A)) . Gjj is the gauge transformation associated 
with U such that GjjU leaves {wy} invariant. 

There is an important subgroup of PSG, Invari- 
ant Gauge Group (IGG), which is composed of all 



the pure gauge transformations in PSG: IGG = 
{{Wi}\W,u,jW] = u,j,W, e SU{2)}. One can always 
choose a gauge in which the elements in IGG is site- 
independent. In this gauge, IGG can be global Z2 trans- 
formations: {Wi = r", Wi = — T^}, the global C/(l) trans- 
formations: {W, = e^^^\e e [0, 27r]}, or the global SU{2) 
transformations: {Wi = e'^"'^, 6* e [0,27r],n e 5'^}, and 
we dub them Z2, U{\) and SU{2) state respectively. 

The importance of IGG is that it controls the low- 
energy gauge fluctuations. Beyond mean-field level, fluc- 



tuations of Uij and On need to be considered and the 



mean-fleld state may or may not be stable. The low- 
energy effective theory is described by fermionic spinon 
band structure coupled with a dynamical gauge field of 
IGG. For example, Z2 state with gapped spinon disper- 
sion can be a stable phase because the low-energy Z2 dy- 
namical gauge field can be in the deconfined phase2ii2^. 
But for a U{1) state with gapped spinon dispersion, the 
C/(l) gauge fluctuations would generally drive the system 
into confinement due to monopole proliferation^^, and 
the mean-field state would be unstable. And an SU{2) 
state with gapped spinon dispersion should also be in 
the confined phase because there is no known IR stable 
fixed point of pure SU{2) gauge theory in 2+1 dimen- 
sion. Because the purpose of this paper is to search for 
stable spin liquid phases that has a Schwinger fermion 
mean-field description, we will focus on Z2 states. 

li GuU e PSG and g e IGG, by definition we 
have gGjjU e PSG. This means that the mapping 
h : PSG -^ SG : f{GuU) = [/ is a many-to-one map- 
ping. In fact it is easy to show that mapping h induces 
group homomorphism^^: 



PSG /IGG = SG. 



(7) 



Mathematically PSG is an extension of SG by IGG. 

Our definition of PSG requires a mean-field ansatz 
{uij}. With Eq.©, one can define algebraic-PSG which 
does not require ansatz {uij}. An algebraic-PSG is sim- 
ply defined as a group satisfying Eq.©. Obviously a 
PSG (realizable by an ansatz) must be an algebraic- 
PSG, but the reverse may not be true, because sometimes 
an algebraic-PSG cannot be realized by any mean-field 
ansatz. 

To classifying all possible Z2 Schwinger-fermion mean- 
field states, we need to find all possible PSG group ex- 
tensions of the SG with a Z2 IGG. Here SG is the direct 
product of the space group of honeycomb lattice and the 
time-reversal Z2 group. In appendix |X] we show the gen- 
eral constraints that must be satisfied for such a group ex- 
tension. In appendix [Bl using these constraints, we find 
there are in total 160 Z2 algebraic-PSGs on honeycomb 
lattice. And at most 128 PSGs of them can be realized 
by an ansatz {uij}. These 128 PSGs are the complete 
classification of Z2 spin liquids on a honeycomb lattice. 



III. CLASSIFICATION OF Z2 STATES AROUND 
THE U-RVB STATE 

Can one further identify the candidate states for the 
spin liquid discovered in the numerical studyi? The an- 
swer is yes. Numerically the spin liquid phase is found 
close to the Mott transition and it seems to be connected 
to the semimetal phase by a continuous phase transition. 
What are the Z2 Schwinger-fermion states in the neigh- 
borhood of the semi-metal phase? 

Are there Schwinger-fermion mean-field states that can 
be connected to the semi-metal phase via a continuous 
phase transition? This question was firstly discussed by 
Hermele in Ref. [23. Using slave-rotor formalism, it was 
shown that the semi-metal phase can go through a con- 
tinuous phase transition into an SU{2) u-RVB state (also 
termed as algebraic spin liquid (ASL) in Ref. 1271 ) at the 
mean-field level. This SU{2) u-RVB ansatz, in terms of 

/-spinon, can be written as Hmf = tJ2<ij> fLfja^ ^ 
is real and summation is over all nearest neighbor bond. 
The single-spinon dispersion of u-RVB state is similar to 
the electronic dispersion in the semi-metal phase, which 
is composed of four two-component Dirac cones at the 
corner of Brillouin Zone, two from spin and two from 
valley. Physically it is easy to understand u-RVB state 
connecting with the semi-metal phase: At the Mott tran- 
sition, only the charge fluctuation becomes fully gapped 
and the spinon dispersion still remember the semi-metal 
band structure. 

The u-RVB ansatz can be simply expressed as a 
graphene-like nearest neighbor hopping of /-fermions: 
FigUl 



rruRVB _ ,, 



z^ 

<ij> 



JiaJjo!^ 



(8) 



where x is real. Beyond mean-fleld level, the low- 
energy effective theory of u-RVB state is described by 
Nf = 2 two-component Dirac spinous {SU{2) gauge 
doublet) coupled with a dynamical SU{2) gauge field^l, 
i.e. QCD3. In the large- A^/ limit QCD3 has a stable IR 
fixed point with gapless excitations and can be a stable 
ASL phaso2^. When Nf = the pure gauge QCD3 is in a 
confined phas e^^'^° . This indicates a critical Nc and when 
Nf < Nc confinement occura^^. Although no controlled 
estimate of Nc is available, a self-consistent solution of 
the Schwinger-Dyson equations^S, suggests A^^ ~ t- ■ We 
will assume that Nc > 2 and therefore u-RVB state is 
not a stable phase. 

Due to the lack of the knowledge of the confinement 
mechanism, it is difficult to reliably predict the ultimate 
fate of the u-RVB state (or ASL). But one possibility is 
that the strong gauge interaction induces Higgs conden- 
sation which breaks the SU{2) gauge symmetry down to 
Z2 , so that the renormalization group flows into a stable 
flxed point of Z2 gauge theory. Based on this assump- 
tion, we can propose a scenario of a continuous phase 
transition from the semi-metal phase into a Z2 spin liq- 




Semi-metal Mott-transition 



ASL 



FIG. 2: Schematic RG flow of the Mott transition. A repre- 
sents a relevant perturbation of the ASL fixed point, which 
eventually drive the RG flow into a stable Z^ spin liquid fixed 
point. 



uid phase: the critical point is still described by the slave- 
rotor critical theory discussed in Ref.[23. But on the Mott 
insulator side, a dangerously irrelevant operator (for ex- 
ample, can be a four-fermion interaction term) becomes 
relevant and finally drive the RG flow away from the ASL 
fixed point and flow into a stable fixed point of a Z^ phase 
by Higgs mechanism. Here we assume the ASL still de- 
scribes an unstable fixed point with relevant directions. 
This scenario is schematically shown in Figl2l which has 
the same spirit of the deconfined quantum criticality^. 

If this scenario is correct, the mean- field ansatz of the 
Zi spin liquid should be connected to the u-RVB ansatz 
by a continuous Higgs condensation, which breaks the 
SU{2) IGG down to Z2. During this transition, the u- 
RVB ansatz {uff^^} -^ {</^^ + Su^J} and the fey 
amplitudes play the role of the Higgs boson. We define 
a Z2 state to be around (or in the neighborhood of) the 
u-RVB when the Z2 state can be obtained by an infinites- 



,uRVB 



}^{ 



,uRVB 
Hj 



,uRVB 



+ SUij}. 



+ Suij} must be a subgroup of the 



imal change {' 

ThePSGof{' 
PSG of the u-RVB state Eq.®. In appendix [C] we clas- 
sify aU these possible PSG subgroups with the Z2 IGG, 
which allows us to construct all possible Z2 states around 
the u-RVB state. This technique was firstly developed by 
Wen2^. We find 24 gauge incquivalcnt Z2 PSGs as listed 
in Table U in appendix [Cl 

Can these 24 Z2 SL states have a full energy gap? We 
find not all of them can have a gapped spinon spectrum. 
This can be understood starting from a Dirac disper- 
sion of the u-RVB state. To gap out the Dirac nodes, 
at least one mass term in the low-energy effective the- 
ory of a given Z2 state must be allowed by symmetry. 
In appendix [E] we show that only 4 of the 24 Z2 states 
allow mass term in the low energy theory. Thus only 
these 4 states are fully gapped Z2 spin liquids around u- 
RVB state. The other 20 states have symmetry protected 
gapless spinon dispersions. 

These four states are state #16, #17, #19, and #22 in 
Table U in appendix [C] We can generate their mean- field 
ansatzs by these PSGs. We find that up to the 3rd neigh- 
bor mean-field amplitudes U(a,0.-y) as shown in FiglU only 
one of these four states can be realized, which is state 
#f9. As shown in appendix IE 21 mean-field ansatzs up 
to the 3rd neighbor of the other three states actually have 



^^ s-wave 
x^pairing 




0.0 0.2 0.4 0.46 0.6 0.8 1.0 

J2/iJl + J2) 

FIG. 3: Mean-field phase diagram of Ji — J2 model by 
Schwinger-fermion approach. 



a C/(f ) IGG. Only after introducing longer-range mean- 
field bonds can these three states have a Z2 IGG. In par- 
ticular, state #16 requires 5th neighbor, state #17 re- 
quires 4th neighbor and state #22 requires 9th neighbor 
amplitudes, while state #19 only requires 2nd neighbor 
amplitudes. Because the t/U expansion of the Hubbard 
model give a rather short-ranged spin interaction for the 
SL phase found in numerics^ {t/U ~ 1/4), the other three 
states are unlikely to be realized in a Hubbard model on 
honeycomb lattice. 

After choosing a proper gauge, the mean-field ansatz of 
#19 can be expressed as a sublattice dependent pairing 
of the /-spinons, as shown in FigHJ 



Hmf -X E /-/^" + ^^''' E '-t^flJ: 



t 






<<ij»eA 

E e„/j/l4+h-c- (9) 



«ij»eB 



and we term it as sublattice pairing state (SPS). Note 
that 6* 7^ 0, ±7r/2, tt, because otherwise the ansatz has 
U{1) IGG. We propose SPS to be the SL phase found in 



numerics. 



IV. SCHWINGER-FERMION MEAN-FIELD 

STUDY OF THE Ji J2 MODEL ON 

HONEYCOMB LATTICE 

Can SPS be realized in the Hubbard model when 
t/U ~ 1/4, where numerics shows a gapped SL phase? In 
particular, by the Mott transition theory of Hermele^l, 
the u-RVB (or ASL) state is in the neighborhood of the 
Mott transition. Can SPS be more favorable than the 
ASL state? To address this question, we use t/U expan- 
sion of the Hubbard model^ to obtain an effective Ji — J2 
spin model on honeycomb lattice: 



H ^ Ji y^ Si ■ Sj + J2 y^ Si ■ Sj 



(10) 



<u> 



<<ij>> 



where Ji and J2 are the 1st neighbor and 2nd neighbor 
antiferromagnetic couphng. Following Ref. |32| . we find 
up to f^/U^ order, the effective Ji and J2 are: 



Ji = Af/U - 16t^/U^ 



J2 = Atyu^ 



(11) 



Naively plugging in t/U ^ 1/4 gives J2/J1 ^ 1/12. 

We use the variationally mean-field ansatz Eql9l Note 
that this mean-field study is biased towards spin disor- 
dered ground state. For example, we do not include Neel 
order which is known to be the ground state at J2 =0, 
and we also do not include the spiral spin order which is 
found by semiclassical study of Ji — J2 model^i^i. The 
purpose of the current mean-field study is to understand 
whether a gapped spin liquid can be more favorable com- 
pared to the gapless ASL state when J2 is tuned up and 
frustration becomes important. 

By minimizing the mean- field energy in Eq.Q, the 
phase diagram of Ji — J2 model is obtained and shown 
in FiglH where we fix Ji + J2 = 1 and Emf is sealed 
from Eq.© by 8/3. We find that when J2/J1 < 0.85 
(or J2/(Ji -I- J2) < 0.46), the ground state is u-RVB(or 
ASL) state: x 7^ and A = 0. When J2/J1 > 0.85, 
the ground state is an s-wave pairing state: Xi ^ 7^ 
and 0^0. The s-wave pairing state opens an energy 
gap for spinous but has remaining U{1) gapless gauge 
fluctuation. Due to monopole proliferation^^ the s-wave 
pairing state is not a stable phase. In this mean-field 
study, the gauge fluctuations are not considered and this 
is the reason why we find s-wave pairing state as a ground 
state. Taking gauge fluctuations into account, the likely 
fate of the s-wave pairing state is that 9 becomes nonzero 
and the Z2 SPS state is realized. 

We propose to study the Ji — J2 model by Gutzwiller 
projected wavefunction variational approach^ because it 
can be viewed as a method to include the gauge fluctu- 
ation. We leave this projected wavefunction study as a 
direction of future research, which may realize SPS as the 
ground state. Projected wavefunctions are also classified 
by PSG, so the present work also provide guideline for 
the search of ground states in the projected wavefunction 
space. 



V. DISCUSSION 

In this work we completely classified the Z2 mean- 
field states in the Schwinger-fermion approach. Using 
physical argument, we identify a single state: SPS, as 
the possible spin liquid phase found in the recent Quan- 
tum Monte Carlo study of the Hubbard model on a 
honeycomb latticei. SPS is in the neighborhood of the 
scmimetal phase and we propose a scenario for the con- 
tinuous transition connecting the two phases. 

In our mean- field study of the J1 — J2 model, the s-wave 
pairing state is realized for a fairly large J2, correspond- 
ing to a fairly large t/U ^ 0.44. A higher order spin-spin 
effective interaction such as the 6-spin ring exchange term 



and/or a more careful projected wavefunction study may 
realize SPS phase for a smaller t/U. 

In a recent work^^, Wang study the Z2 mean-field 
states in the Schwinger-boson approach, and identify a 
zero-flux SL state, which is naturally connected to a Neel 
ordered state by a potentially continuous phase transi- 
tion. Whether the SPS found in the present work is 
related to Wang's result is unclear. And we leave the 
possible continuous transition from SPS to the Neel or- 
dered phase as a subject of future research. 

YR thanks Ashvin Vishwanath and Fa Wang for help- 
ful discussions. YML thanks Prof. Ziqiang Wang for 
support during this work under DOE Grant DE-FG02- 
99ER45747. YR is supported by the start-up fund at 
Boston College. 



Appendix A: General conditions on projective 
symmetry groups on a honeycomb lattice 

As mentioned in section |IT1 SG on a honeycomb lattice 
is generated by time reversal transformation T, trans- 
lations along 0,1,0.2'- Ti , T2 , plaquctte-centcred 60° Cg 
rotation, and a horizontal mirror reflection cr as shown 
in Fig[TJ In the present problem, the symmetry group 
SG can be represented as 



SG^{U 



rpUT , rp^Ti ^ rp'^T2 f^'^Cg 



r-^'^} 



where vti , 1^X2 ^ ^ ^-nd vt^Vo- G .Z2, vce ^ ■^e- since the 
generators satisfy 



{C,f 



1 



(Al) 



Here 1 stands for the identity element of SG. To com- 
pletely determine the multiplication rule of this group, 
we need to identify the multiplication rule of two differ- 
ent generators in an order different from T"^^ • Tj^ ^^ • T2 ^^ • 



UT = TU {U = Ti,T2,C6,(t) 


(A2) 


T1T2 = T2T1 


(A3) 


CqTi = T2CQ 


(A4) 


CeT2 = T,-'T2Ge 


(A5) 


ctTi = Ticr 


(A6) 


(TT2 = TiT^^o- 


(A7) 


(tCq = Cq (t 


(A8) 


The above relations can be written in an 


alternative 


way 





T^ = cr2 = {Cef = 1 (A9) 

TC/T"V-i = 1 (t/ = Ti,r2,C6,cr) (AlO) 

(All) 
(A12) 
(A13) 
(A14) 
(A15) 
(A16) 



TiTiTj^^T^^ = 1 
T2 CeTiCg = 1 

T^^a-Ticr-^ = 1 

crCgcrGQ = 1 



which determines the inverse of all the group elements. or more specifically 

As introduced in section|lTl the mean- field ansatz {uij} 
of a spin liquid is invariant under the action of any ele- 
ment GjjU of a projective symmetry group (PSG). The 
multiplication rule of the symmetry group would imme- 
diately enforce the following constraints on a PSG by its 
definition: if U1U2 = U3 then 



[Gu^{UiU2{i))Gu,{U2it))]u,,[Gu,{UiU2{t))Gu,{U2{i))]^ 
= [Gu3{U3{t))]u,j[Gu3{U3{i))]\ ViJ (A17) 

On the other hand, we know those pure gauge trans- 
formations, under which the mean- field ansatz {uij} is 
invariant, constitute a subgroup of PSG, coined the in- 
variant gauge group (IGG): 

IGG^{W,\W,u,,w] =u,j, W,eSUi2)} (A18) 



Gt'{ 



Therefore from (jA17p we have the following constraints 
on the elements of a PSG 

[Gu,u,{UiU2{i))yGu,{UiU2{t))Gu,{U2{i)) = G e IGG 

The above condition holds for any two group el- 
ements Ui,U2 of SG. Similar with SG, we can 
choose a set of generators in any given PSG: 
{Gt,Ti,Gt,T2,GtT, Gc.Ce, G^ct}. Any given element 
in PSG can be written in the standard form: 



GuU ={GtTY^ ■ {Gr^Tiy^i ■ {Gt,T2Y^^ 



(A19) 



Since the multiplication rule of SG on a honeycomb lat- 
tice is completely determined by (|Aip - (jA8p . or equiva- 
lently (|A9[) - (|A16|) . the only independent constraints on 
the PSG generators are the following: 



[Gr(J)] e/GG, (A21) 

G,(ff(z))G,We/GG, 

Gc,{C^Ht))Gc,{C^^ii))Gc,{Gl{i)) 

■Gc,{Cl{i))Gc,{G^{i))GcS) G IGG, 

'T^^T,{{))G^I{T,{i))Gt,{T,{^)GtS) e IGG, 

G^^{Ti{i))G^'{Ti{i))GT,{Ti{i))GT{t) G IGG, 

G^^{T2{i))Gj,'{T2{i))GTAT2{i))GT{i) e IGG, 

G^^{T2ii))Gc,iT2ii))GTATiG^\^))Gcli^) e IGG, 

G^I{T,{{))Gc,{T,{i))Gt,{C^^T^{i)) 

■G^l{C^\i))G-cl{i)€lGG, 

G^l(G6-l(^))Gcl(^)GT(^)Gc«(^) £ IGG, 

G^I{T^{i))G„{T^{i))Gt,{Tict-\{))G„\i) e IGG, 

G^^{T2{t))G^{T2{t))GT,{aT2{t))G^^{(T{t))G^\t) e IGG, 

G^(^)Gce(<T(^))G^(^C6(^))Gca(G6(^)) G IGG, 

Gr'(^(j))G^H*)GT(*)G^(i) e IGG. 

Above are all the general consistent conditions to be sat- 
isfied by the generators of a PSG on a honeycomb lattice. 
We will use {xi, X2, s) to label a site i in a honeycomb 
lattice, where xi,a;2 are the coordinates of the unit cell 
in basis 0,1,0,2 and s = 0, 1 for A and B sublattice re- 
spectively. For convenience, we summarize the coordi- 
nate transformation of all the generators in the symmetry 
group on a honeycomb lattice as follows: 



Gfi 



T : {xi,X2,s) -> {xi,X2,s), 
Ti : {xi,X2,s) -)■ {xi + l,X2,s), 
T2 : {xi,X2,s) ^ {xi,X2 + l,s), 
: {xi,X2,s) ^ {Xi+X2,-X2,l- s), 
: (a;i,a;2,0) -> (1 -X2,a;i +yi - 1,1) 
(xi,X2,l) -^ (-2-2, a;i +yi,0) 



(A22) 



[GtTY e IGG (A20) 

(G^o-)2 e IGG 
(GcCeY € IGG 

{Gt,TiY\Gt,T2Y\Gt,Ti){Gt,T2) £ IGG 
{Gt,Ti)-\GtTY\Gt,Ti){GtT) e IGG 
{Gt,T2Y\GtTY\Gt,T2){GtT) e IGG 

{GT,T2r\Gc,GQ){GT,T^){Gc,G^r^ e IGG 

{GT,TiY^iGc,Ge){GT,Ti){GT,T2)-'{Gc,CeY' e IGG 

iGTTYHGceCeYHGTTYGcGe) € IGG 

{GT,T^r\G^cT){GT,T^){G^cTY^ G IGG 

{GT,T2YHG.cT){GT,T{){GT,T2r\G^aY^ G IGG 

(G^o-)(GceG6)(G^rT)(Gc«G6) e IGG 

iGTTYHG.TYHGTT){G^(T) G IGG 



Appendix B: Classification of Z2 projective 
symmetry groups on a honeycomb lattice 

As discussed in sectionllll the problem of classifying all 
possible Z2 Schwinger-fermion mean-field states is math- 
ematically reduced to finding all possible PSGs. Let us 
firstly find all algebraic PSGs. 



1. General discussions 

In the case of Z2 spin liquids, the IGG of the cor- 
responding PSG is a Z2 group: IGG = {±t°}. The 



constraints listed in appendix |X] now becomes 

[Gri^f^VTr'', (Bl) 

G^{fTii))G^ii)^f^^T°, (B2) 

Gce{C^'{^))Gc,{C,\i))Gc,{C|{^)) (B3) 

■Gc,{Gi{{))Gc,{Ce{i))GcA^) = VCst", (B4) 
G^^{T,:'T,i^))G^^{T,{i)) 
•GTi(Ti(i))GT,(i) = 7712^0, (B5) 

G^^{T,{i))Gj,'{Ti{i))GT,iTi{i))GT{i) = r/iTrO(B6) 

GTl{Ui))Gc.{T^{i)) 
■Gt, {T,G^'{i))Gcl{i) = ?7c«ir", (B8) 

G^^i (Ti (z)) Gc« (Ti (i)) Gt, (Gg-^Ti (*)) 

■G^l{G^'{i))Gcl{i)^Vc,2T^, (B9) 

G^i(G6-i(z))G^i(z)GT(i)Gce(*) =r;c,Tr", (BIO) 
G^^i(ri(i))G.(ri(*)) 
■Gt, {T,a~\i))G^\i) - r/^irO, (Bll) 

G^l{T^{^))G^{T^{^))GT,{cTT^{^)) 

•G^^i(<T(^))G,i(z) = r;^2rO, (B12) 

G.(z)Gce(T(z)) 
•G^(<TG6(i))Gc,(G6(z)) = V.cer\ (B13) 

Gyi(<T(i))G-i(i)GT(«)G^(z) = t^^ttO. (B14) 

where all the 77's take value of ±1. Not all of these condi- 
tions are gauge independent. Because we can re-choose 
the gauge part of generators such as Gti , Gt2 • ■ • by mul- 
tiplying them by —t^ (an element of IGG), only those 
conditions in which the same generator shows up twice 
are gauge independent. Wc can use this gauge depen- 
dence to simplify these conditions. Because GxiiGr^) 
only show up once in the equation of r]Cei{riCe2) , we can 
always choose a gauge such that jycgi = VC62 ~ 1- AH 
other 77's are gauge independent. 

In the following wc will determine all the possi- 
ble PSG's with different (gauge incquivalent) elements 
{Gu{i)}- These different PSG's characterize all the dif- 
ferent type of Z2 spin liquids on a honeycomb lattice, 
which might be constructed from mean- field ansatz {uij}. 

First notice that under a local SU{2) gauge transfor- 
mation Uij —J- WiUijW^ , the PSG elements transform as 
Guii) — > WiGuii)W^_i,.y Making use of such a degree 
of freedom, we can always choose proper gauge so that 



Gti{xi,X2,s) = GT2iO,X2,s) 



a;i,.T2 e Z. 



Now taking (jBSP into account, we have Gt2{{xi + 
l,a;2,s}) = 77i2GT2({a;i,X2, s}) and therefore 



GTiixi,X2,s) 
Gt2{xi,X2,s) = 



(B15) 



Meanwhile, from (|B1[) . (|B6I) and (|B7|) we can immedi- 
ately see that tjit = 772T = 1, and the gauge incquivalent 
choices of Gxii) are the following 

GTixi,X2,s)=gT{s)= I '* ' (B16) 

where rjt = ±1. 

As discussed earlier, we can always choose a proper 
gauge so that rycei = ^C6 2 = 1- Then from conditions 
and (IB9I) we see that 



GcA^Y,X2,s)=T]^^ " gce{s) (B17) 

similarly from conditions (JB11|) and ()B12p we have 

G.(xi,X2, s) = v7iV72Vit''''^^'9A^) (B18) 



where gceis),g^{s) e SU{2). Note that (1B2| an d (|B13jl 
give further constraints to the above expression (jB18[) : 

Vo-i = Vo-2 = ?7i2 (B19) 

Now we see the elements of PSG can be expressed as 

GtAxi,X2,s) = t" (B20) 

Gt2{xi,X2,s) =7712'^° 

GT{xi,X2,s)=gT{s) (B21) 

ri I \ 2;i2;2+a:i(a:i — 1)/2 / n 

Gc6(a;i,a:^2,s) = '7i2 9cAa) 

G„{xi,X2,s) = r]i2~^ "" yrris) 



Consistent conditions jlSI, (|B4|, (|BT0| . (|BT3| and (|BT41) 
correspond to the following constraints on SU(2) matri- 
ces gceis), 9ct{s): 



[9Ce{s)9Ce{'^ - s)f = VCeVl2r°, 
9T{s)9Ce{s) = 9Ce{s)9T{l ~ s)lJCeT 

9T{s)g^{s) = g^{s)gT{l - s)t]^t 



(B22) 



Ac„T°, V^Ce = 1 



^ iris ■ T, 7]^Ce 

where Acg = ±1 and fis is a unit vector. 



2. A summary of 160 different PSG's 

Below we summarize all the 160 possible PSG's 
obtained through solving (JB22|. We use capital Roman 
numerals (I) and (II) to label gx = yfr'^ and gx = ir^ 
respectively. Roman numerals (i) and (ii) are used to 
label 77C6T = ±1 respectively. (A) and (B) are used to 
label 77ctC6 = il respectively. Finally (a) and (/3) are 
used to label ?7ctt respectively. 



„a:i 
'/l2~ 



(I) 9o 



Vlr": 



It's easy to see that tjcqT ~ VrrT = Vt from (jB22p . so 
there is the only possibility among (i) and (ii). 
(A) g^is) ^ X-^^gclil ~ s): 
we have Xce = il(TilCeVi2 and 



5<t(1) ^ VTVCeVl2T° ■ 



(B23) 



This represents 2'^ = 16 different PSG's in the class (I) (A) 
since ?7t, r/Ce, '7<t,??12 = ±1- 

(B) g^{s)gce{i - s) = ihs ■ f: 

Choosing a proper gauge (so that gc'ei^) = '''") ^^ have 



(B24) 






where V'3 = 0,±27r/3 stand for the multiples of 2tt/3 
mod 2n. There are 2" x 3 = 48 different PSG's in this 
class (I)(B). 



(II) 5T(s) = iT3: 

(i) ??C6r = 1: 

(A) g^s) = X^CedcH^ - s)- 

in this case Ace = Vo-VCeVi^, so we have 

.9c«(0)=t0, 

ff<T(0) =5C6(1) = VCeVl2T°, 
5<t(1) ^ VcTVCeVl2r° ■ 



(B25) 



there are 2^ = 8 different PSG's in the class (II)(i)(A). 
(B) g^{s)gce{^ - s) = iris ■ f: 
(a) ija-T = 1, «-e. [go-is), T^] = 0: 
here we have 



go-(O) = -ir/o-r^, 
3<x(l) = ir3. 



(B26) 



here we must have r]crT = ^1, ^Ce = 'yo-'yce ^12 and 

gcM - iT\ (B28) 

506(1) = -iVCem2T\ 
ffo-(O) = iVCeVl2T^, 



ff<x(l) 



-l??<TfyC6'7l2T 



there are 2^ = 8 different PSG's in the class (II)(ii)(A). 

(B) go-(s),9c6(l - s) = i"s • ■^^ 
(a) r?o-T = 1, «-e. [go-(s),T3] ^ q. 
here we have 



50^(0) = iT^, 
5C6(1) = -iVCeVi2T^e''f'^'''' 



(B29) 



there are 2^ 
(II)(i)(B)(«). 

(/?) VcrT = -l,t-e.{gA^s),T^} = 
here we have 



5.(0) = rO, 
.9.(1) -^.r". 

there are 2^ x 3 = 24 different PSG's in the class 
(II)(ii)(B)(a) since ^3 = 0,±27r/3. 
(/3) 77.T = -1, J-e. {ffcr(s),T3} = 0: 
here we have 

.9Ce(0) = iri, (B30) 

ffce(l) = -ir/c.m2Tie"^^-', 

.9.(0) = iTi, 

5.(1) = -iv^^- 

there are 2^ x 3 = 24 different PSG's in the class 

(II)(ii)(B)(/3) since V3 = 0,±27r/3. 

To summarize, above are the 160 different (algebraic) 
PSG's with IGG = {±t"} on a honeycomb lattice. They 
represent different Z2 spin liquid states on a honeycomb 
lattice, which possess all the symmetries of the honey- 
comb lattice generated by {T,Ti,T2,cr, Ce}. We also 
want to emphasize that any solution to the set of equa- 
tion (jBip - (jB14p may look different, but it will be gauge 
equivalent to one of these 160 PSG's. 

On the other hand, such a (algebraic) PSG really cor- 
responds to a spin liquid if and only if it can be realized 
by a mean- field ansatz {uij} on a honeycomb lattice^^. In 
fact, not all of these algebraic PSGs can be realized by an 
ansatz. After the time-reveral transformation, the mean 
field amplitude changes signS^: T(u,y) = —Uij. Gauge 



8 dirterent PSG s m the class transformation Gt must change the sign again 



-Uij = GTii)uijGTU)^ 



(B31) 



If in an algebraic PSG, Gxii) = t'^ independent of site, 



(B27) 



5C6(0) = tO, 

gCei^) = VCeVl2e''^'^\ 

5.(0) = -i?7.T\ 
5.(1) = ir\ 

there are 2^ x 3 = 24 different PSG's in the class 

(II)(i)(B)(/3) since Va = 0,±27r/3. 
(ii) VCeT = -1: 
(A) 5^(s) = A^^gci(l-s): 



Uij must vanish. 



Clearly at least 32 algebraic PSG's among the total 
160 types cannot be realized by any mean- field ansatz 
{uij}. These are the PSG's with Gxii) = gxis) = t° 
in the class (I)(I)(A)&(B). Since under time reversion T 



ij 



GTii)uijGl{j) 



we requue — m 

the vanishing of all bonds {i 



H], 



this leads to 
0} in the mean-field 
ansatz. Therefore, there are at the most 128 possible Z2 
spin liquids that can be realized by a mean-field ansatz 
on a honeycomb lattice. 



Appendix C: Classification of Z2 projective 
symmetry groups around u-RVB ansatz 

In this section we focus on those Z2 spin hquids near 
the u-RVB state, which is discussed in section Hill These 
Zi spin hquids are plausibly connected to a scmimetal 
through a continuous phase transition. The u-RVB state 
is realized by the following ansatz: 



(-!)«• ixr° 



(CI) 



its mean-field bond is only nonzero between nearest 
neighbors < ij >, which have different sublattice indices 
;. By definition, its PSG has the following 



S^ = 1 

form: 



Gti{xi,X2,s) = gi, (C2) 

Gt2{xi,X2,s) = g2, 

Gt{xi,X2,s) == (-1)*5T, 

Gce{xi,X2,s) = {-lygcs, 

G(^{xi,X2,s) = {-lyga-. 

where gi,g2,gT,gce.i9cr € SU{2). To find out those 
Z2 spin liquids around such a u-RVB state, wc need to 
trace those PSG's with IGG = {±t°} that looks like 
(jC2[) . Consistent conditions (|Bll) - (|B14p now corresponds 
to constraints on the SU(2) matrices {gi, 52, .9t, gce 1 9o-}'- 

c ^0 

= 42TT , 

C ^0 
= ?<t2T , 



9i^92^9i92 


= ei2r°, 


2 

9t 


9i^9t^9i9t 


= eiTr°, 


92^9t^929t 


92^9Ce9l9cl 


= eceir°. 


9T^9Ce9l92^9de 


9T^9Ce9T9Ce 


= ^CeTT , 


9Ce 


9i^9t9i9^^ 


= e^ir°. 


92^9ct9i92^9ct^ 


9cT9Ce9T9C6 


= Co-Ce'T , 


9t^9^^9t9o- 



^.tt" 



C<.r°. 



(C3) 



where all ^'s take value of ±1. Again, as discussed in 
appendix [B] we can always make ^^gi = ^C6 2 = 1 by 
choosing a proper gauge. After solving eqs. (|C3p . we find 
out there are 24 gauge inequivalent solutions in total, as 
summarized in TABLE ID In other words, there are 24 
different Z2 spin liquid around the u-RVB state, suggest- 
ing that they are promising candidates of the spin liquid 
connected to a scmimetal on honeycomb lattice through 
a continuous phase transition. 



# 


qt 


3<T 


gce 


31 


ff2 


1 


r" 


r" 


r" 


r^' 


r" 


2 


r" 


r" 


ir^ 


r" 


r" 


3 


t" 


r" 


iT^ 


gi2,r/3r^ 


g-i27r/3Ti 


4 


r" 


ir^ 


ir^ 


r" 


r" 


5 


r" 


ir^ 


ir^ 


r" 


r" 


6 


r" 


ir^ 


ir^ 


r^' 


r" 


7 


r« 


ir^ 


gi./6r^ 


r» 


r° 


8 


t" 


ir^ 


ei-/3r^ 


r" 


rO 


9 


t" 


ir^ 


ir^ 


gi27r/3r-* 


g-i27r/3T^ 


10 


r« 


ir^ 


gi2,r/3rl 


i(^-V'F) 


n^f "75 "73'' 


11 


ir' 


r" 


r" 


r" 


r" 


12 


ir' 


r" 


ir^ 


r" 


r" 


13 


ir-^ 


r" 


w' 


r» 


r" 


14 


ir^ 


r" 


■it' 


gi2^/3r-' 


g-i27r/3T-' 


15 


ir' 


ir^ 


r" 


r" 


r" 


16 


ir' 


ir^ 


ir^ 


r" 


r" 


17 


W' 


ir^ 


ir^ 


r" 


r" 


18 


ir^ 


ir^ 


W' 


gi2,r/3r-' 


g-i27r/3T-' 


19 


ir^ 


ir^ 


W' 


r" 


r" 


20 


W' 


ir^ 


\t' 


r" 


r" 


21 


W' 


ir^ 


r" 


r" 


r" 


22 


ir' 


ir^ 


ir^ 


r" 


r" 


23 


ir^ 


it' 


giT/6r^ 


r" 


r" 


24 


ir^ 


it' 


giT/3r-' 


r" 


r" 



TABLE I: A summary of all 24 different PSG's with IGG = 
{±r } around the u-RVB ansatz. They correspond to 24 
different Z2 spin liquids near the u-RVB state. 



elements must transform mean-field bond Uij into itself 
Uij or its Hcrmitian conjugate uj- = Uji. 

As a special case, the identity element 1 always trans- 
form a bond into itself: correspondingly in PSG the IGG 
elements (e.^. r° for a Z2 ansatz) always transform any 
bond Uij into itself. This is nothing but the definition of 
invariant gauge group (IGG). 

Now we need to look at nontrivial symmetry group 
elements which transform two lattice sites (connected by 
the bond) into itself or into each other. Without loss of 
generality, we consider the following bond 



\X\^X2-,s) = U(xi,a;2,s)(0,0,0) 



1. Regarding time reversal T 



(Dl) 



Appendix D: Consistent conditions on the 
mean-field ansatz \uij\ on a honeycomb lattice 

In this section wc derive the consistent conditions on 
an arbitrary mean-field bond u,,,, which realizes a spin 
liquid with a certain PSG on a honeycomb lattice. The 
basic idea is to find all possible symmetry group elements 
that transform the two lattice sites {i, j} into itself {i, j} 
or into each other {j, i}, so that the corresponding PSG 



Any clement of the symmetry group can be written as 



7-7- _ rpVT rnl^T, m^'T^, _j/„ r^^Ce, 



-'l ' ^2 



■ c: 



(D2) 



First we study the consistent conditions from time re- 
versal transformation T and then turn to other group 
elements. 

Notice that time reversal T doesn't change anything 
except the sign of bond: 



Gt{^u,,[Gt{])]^ 



(D3) 
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so this bond must satisfy the foUowing constraint: 

GT(a;i, 2^2, s)(.Ti, 2:2,5) (D4) 

= -(xi,.T2,.s)Gt(0,0,0) 



for V x, xi, 2:2 S Z. 



4. An example: mean-field ansatz {uij} of Z2 spin 
liquids near u-RVB state 



2. Conditions on a bond within the same 
sublattice: s = 

First we study s = case, i.e. a bond within the same 
sublattice. Since both mirror reflection cr and 7r/3 ro- 
tation Cg will change the sublattice index s while the 
translations Ti,T2 don't, we must have an even number 
of reflection and rotation, i.e. v^- + vc^, — mod 2 to 
transform the bond to itself (or its Hermitian conjugate). 

^From (|A22p it's easy to check the 5 nontrivial ele- 
ments consisting of {cr, Cg}: 

C|(xi, 2:2,0) = (1 -xi -X2, 2:1,0), 

Cg"^(xi, 2:2,0) =: (2-2,1 -2:1 -.T2,0), 

crC6(xi, 2:2,0) = (2-1,1 -2:1 -yi,0), 
crC| (2:1, 2:2,0) = (1 -.Ti -2:2,22,0), 

(tC6-^(xi,.T2,0) = (2:2,21,0). (D5) 

In order that the bond goes back after some translations, 
it's straightforward to find out all the consistent condi- 
tions on such a bond: 



T2"^crC6 : (-22, 2, 0) -^ (-22:, 2:, 0) 



(D6) 



T^ 



tC&: (0,2:,0) ^ (0,2:,0)t 



T^^crCl : (2, -22:, 0) -^ (x, -2x, 0) 

T^-^<tCI: (x,0,0) -> (x,0,0)t 

(tCq^ : (x, X, 0) -> (x, X, 0) 

TfT-'^crCfr^ : (x, -x, 0) -^ (x, -x, 0) + 



for V X e Z. 



3. Conditions on a bond connecting different 
sublattices: s = 1 

In the s = 1 case, such a bond connects different sub- 
lattices. So only an even number of reflection and ro- 
tation, i.e. Vcr -|- vcf, — mod 2 might transform the 
bond to itself, while an odd number of reflection and ro- 
tation, i. e. Va- -t- vc^ = 1 mod 2 can transform the bond 
(xi,X2, 1) into its Hermitian conjugate (xi,X2, 1)^. 

It's straightforward to obtain the following conditions 
on the bond (xi,X2,l) = W(xi,x2,i)(o,o,o): 



tC« 



(-2x,x, 1) -^ (-2x,x,l)1' 
(x + l,x, 1) ^- (x + l,x, 1) 



(D7) 



Tf^^-^T^+VCg-^ : (-2x- l,x,l) ^ (-2x- l,x,l)t 



-2x-2rpx+\ ri-'i' 

Tl'-^T^-'Cl ■ (X1,X2,1) ^ (xi,X2,l)t 
T^^crCl : (x, -2x, 1) -^ (x, -2x, 1) 



rpX — \rpX — \ 

T2 ctC^ : 



■Cl: (x + l,x,l)^ (x + l,x,l)t 
(-2x - 1, X, 1) ^ (-2x - 1, X, 1) 



To demonstrate the above consistent conditions, let's 
take a look at how they determine the mean-field ansatz 
\uij\ of any Z2 spin liquid near u-RVB state, with PSG 
generators (|C2p . 

Considering time reversion T we immediately have 



gT(2:i,X2,s) = -(-l)''(xi,X2,s)gT 



(D8) 



In other words, the bond connecting two sites be- 
longing to the same (different) sublattice(s) anti- 
commutes(commutes) with qt. 

For the nearest neighbor (n.n.) bond Uq = (0, 0, 1) we 
have xi = X2 = 0,s = 1. Conditions (|D7|) and (|D8|) 
immediately lead to 






(D9) 



For 2nd n.n. bond u^ = (0,1,0) we have xi = = 
s,X2 = 1 and (|D6|) . (|D8|) lead to 



(DIO) 



{u/3,ffT} = 

g^gCeU/s = u^g^gce 



For 3rd n.n. bond u-y = {I, 0, 1) we have X2 = 0, xi 
3 = 1. Conditions (iDTl) and (iDSl) lead to 



K,5t] = 



(Dll) 



r,0 T, 



,t„3. 



Bcrgd'^i = u-ig^gd 



Constraints on further neighbors: e.g. 4th n.n. 
(0,1,1), 5th n.n. (1,1,0) and 6th n.n. (2,0,0) can be 
similarly obtained. 



Appendix E: A search of gapped spin liquids near 
the u-RVB state 



In appendix [C] we showed that there are at most 24 
Z2 spin liquids around the u-RVB state, which are likely 
to connect with a scmimctal through a continuous phase 
transition. In this section we search for those states with 
spectral gaps among the 24 spin liquid ansatz. In the 
end we find out most of the 24 states are gapless. More 
specifically, they cannot open up a mass gap through 
any perturbation around the u-RVB state, which has two 
graphenelike Dirac cones in the 1st Brillouin zone. It 
turns out that only 4 of them, i.e. #16, #17, #19 and 
#22 in TABLElIl are gapped spin hquids near the u-RVB 
state. 



11 



Symmetry-allowed masses in a graphenelike 
u-RVB state 



We start from the low-energy effective Hamiltonian 
of the u-RVB state, which is described by a mass- 
less 8-component Dirac equation. These 8 compo- 
nents contain 2 spin indices (labeled by Pauli matri- 
ces {t*}), 2 sublattice indices (labeled by Pauli matri- 
ces {/i'}) and 2 valley indices (labeled by Pauli matri- 
ces {i^^})- Just like graphene, the two valleys are lo- 
cated at K and K', i.e. the vertices in the honeycomb- 
shaped 1st Brillouin zone. Following the convention 
shown in FIG. [I] the momentum of these two cones are 
K = ^bi + '^b2 and K' = ^h + ^h respectively, 
where {h = (VS, -l)/\/3a, h = {0,2)/V3a} arc the 
reciprocal lattice vectors corresponding to lattice vectors 
{ai = (a,0), a2 = (l,V3)a/2}. 

Expanding the mean-field Hamiltonian of a u-RVB 
state with u„ = ir*^ (here k = — ^(fc^, k„) = kibi +^262) 



ir" (herek= ^(^^'^; 



VJ 



H,, 



/,t ,/,t 



luRVB = i(V'k,A'V'k,B)- 

-r"(l + e-''== +e'('=i-'=^)y 



^0(^^gifc2 _^gi(fc2-fcl)) 







around K and K' we immediately obtain the Dirac equa- 
tions 










T [ky IfCa-j 



T°{ky + ifc^y 



r"(fc'-ifc- 



V'k,B 

i'k',A 

V'k'.B 



[r^ik'y+ik'J 

Defining the following 8-component spinor: 

*''^(V'J,A,'^k,B>k',B>k',A) (El) 

we can write the above effective Hamiltonian of u-RVB 
state as 

iJ = *V^(/^^9,. +/ii9j,)®r°®z/°* (E2) 

Therefore only those mass terms M ^ jt" 



a mass gap can be opened in the Dirac spectrum. In 
the following we study how the mass term changes under 
the action of symmetry transformation such as spin ro- 
tations, time reversal T and translations Ti,T2 etc. The 
physical symmetry of a spin liquid state realized by mean- 
field ansatz only allow those masses that are invariant un- 
der the corresponding PSG. If a PSG already rules out all 
possible mass terms M ~ ir'®T°-®v'^ , a,b = 0, 1, 2, 3, we 
conclude the corresponding spin liquid realized by mean- 
field ansatz is gapless. 

First we work out the transformation rules of Dirac 
spinor "^ and M under a PSG. We focus on the 24 PSG's 
near the u-RVB state with the form (jC2l) as summarized 
in TABLE n 



a. Spin rotations 

It's straightforward to show that a spin rotation along 
z-axis by 26 angle is realized by 



*^e'''* 



(E3) 



while a spin rotation along y-axis by tt angle is realized 
by 



* -^ ir-^V^* 



(E4) 



Apparently Sz rotations leave the mass term invariant, 
while under Sy rotations by tt the mass term transforms 
in the following way 



1 ^ ,,1 ^ ^2 i\,tT2 ^ ,,1 ^ ,,1 



M -^ -/i^ (^ly'- ® t'^M^ t^ ® ^^ ®v 



(E5) 



Since the mass term is invariant under spin rotations, 
its allowed form as seen from the above constraint can 
only be 



m'^^^ = ^J?■ ® v^ ® t" , a ==1,2, 3 



M^g'^ = 11^ ® v^ ® T° , 6 = 0,1,2 



Time reversal T 



(E6) 
(E7) 



Since a mean-field bond uij becomes 
— {—\Y'gTUijgrp[—\Yi under the time reversal transfor- 
mation in a PSG (|C2p . clearly T is realized by 



* 



M - 






(E8) 



so the mass term is invariant under time reversal T if 



M 



-gr' 



I fi^ (g) v^Mg]. ig) n^ igjiy^ 



(E9) 



10 spin liquids near the u-RVB state, i.e. #1-#10 in 



TABLE U has gr 



In these cases, mass terms 



a, 6 = 0, 1, 2, 3 satisfy that {H, *tM*} = so that M 



(a) 



1,2,3 will violate transformation rule (|E9p . 



and the only allowed masses are M^ 

The other 14 spin liquids around u-RVB state (#11- 
#24 in TABLE m are characterized by gr = ir^. In this 
case the allowed masses are Mj^ ' , M^ ' and M^ ' , M^ ' . 



c. Translations Ti, T2 

Under translations Ti,T2 in a PSG ((£2]) the 8- 
component spinor changes as 



T2 









(ElO) 
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since K • ai^2 = T^ and K' • ai^2 = ^^- ^^ order for 
the mass term to be invariant 



,3 






■)giMe^'—' 



'9l 



>9l 



(Ell) 



the symmetry-allowed masses can only be: 
1,2,3 if 51 =g2 = T 



M^g^ and A/jf ^ if gi 



0. 



.92 



,i2Tr/3T-'. 



M^°^ for the special case #10 in TABLE H 

Combining conditions (jE9p and (|E11[) we can see that 
{M^ , & = 0, 1, 2} are not allowed by symmetry in any 
of the 24 spin liquids near u-RVB state. In the following 



study we will focus on masses M 



(a) 
A ' 



1,2,3. 



Reflection cr 



Similar to time reversal T, under reflection along x- 
axis the spinor transforms as 

* ^ /i^ -g^® /^^(8)i^^* = -ig^® /^^ 1^^* (E12) 

The mass term is invariant under reflection cr if 

M = go- ®[i^® i^^Mgl ® i? ® v^ (E13) 



The symmetry-allowed masses are: 
none if ga 



rO. 



M 



(a) 
A ' 



a ^ & if g^ = iT^ 



e. 7r/3 rotation Ca 
Under Cg, i.e. a rotation by 7r/3 the spinor transforms 



as 



* 



9Ce 



^2 2 ' 



(E14) 



The mass term is invariant under reflection cr if 



J 



.i-g-M 



-v^„l _ l,,^\ 



The symmetry-allowed masses are: 

none if gcg = t", e'^"^ with 9 y^ mod tt/2; 



(E15) 



# 


Ua 


Up 


u., 


w« 


We 


9th n.n. (1,2,0) 


16 


ir'J 


{r'^} 


ir" 


ir" 


{r\r^} 




17 


ir" 


r' 


ir" 


{ir^r^} 






19 


{ir",r'^} 


{r'.r'} 










22 


ir" 


t' 


ir" 


ir^' 


t' 


{-\-n 



TABLE II: Symmetry-allowed mean-field ansatz of the 4 pos- 
sible gapped spin liquids near the u-RVB state. We follow the 
notations for mean-field bonds in appendix [D] We only sum- 
marize the mean-field bonds that are necessary to realize a 
gapped Z-i spin liquid. Ellipsis represents those longer-range 
mean-field bonds unnecessary for a Z2 spin liquid, which are 
not listed in this table. Up to 3rd n.n. mean-field bonds 
\uoL, Up, u.y}, only one Z2 spin liquid, i.e. #19 can be realized 
on a honeycomb lattice. 



u-RVB state. Only the following 4 spin liquids near the 
u-RVB state can obtain an energy gap in the spectrum 
through adding a syinmetry- allowed mass term: 



(1,2) 



(2) 



#16 with two symmetry-allowed masses M\ 

#17 with one symmetry- allowed mass M 

(2) 
#19 with one symmetry- allowed mass M\ ' 



if® 



^•^ ® 



(2) 



#22 with one symmetry-allowed mass Af^ 

/3 ^>, ^2 



v^ ® r" 

In fact, as summarized in TABLE|lTl all these 4 gapped 
spin liquids can be realized by mean-field ansatz {wij}, 
which satisfies consistent conditions from the correspond- 
ing PSG as discussed in appendix |D] In the following 
we describe the mean-field ansatz for these 4 gapped Zi 
spin liquids. In the end only one gapped Z2 spin liquid, 
i.e. #19 can be realized by a mean-field ansatz up to 3rd 
n.n. bonds. 



a. Z2 spin liquid #16: up to 5th n.n. bonds needed 

The mean-field ansatz {uy} for Z2 spin liquid #16 
is summarized in TABLE HIl up to 5th n.n. bonds. The 
corresponding spin liquid has a Z2 gauge structure, if and 
only if [up,Ue\ ^ 0, so that the IGG of this mean-field 
ansatz is a Z2 group {±r°}. 

It's straightforward to check that 2nd n.n. bond up = 
/3it1 -f P2t'^ open up a mass gap M ^ ^i^ ®v^ ® {Pit^ + 



P2T')^PlM^^'^+p2M^ 



r(2) 



2. Realizing the 4 gapped Z2 spin liquids near the 
u-RVB state 



Among all 24 spin liquids near the u-RVB states, it 
turns out that there are no symmetry-allowed masses for 
20 of them. In other words, these 20 spin liquids cannot 
open up a mass gap through a perturbation around the 



b. Z2 spin liquid #17: up to 4th n.n. bonds needed 



The mean-field ansatz {uy} for Z2 spin liquid #17 
is summarized in TABLE HIl up to 4th n.n. bonds. The 
corresponding spin liquid has a Z2 gauge structure, if and 
only if [m/3,U5] 7^ 0, so that the IGG of this mean-field 
ansatz is a Z2 group {±r°}. 
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It's straightforward to check that 2nd n.n. bond up 



Pt"^ open up a mass gap M ^ fSfi^ ® v'^ ® t^ 



PMf. 



c. Z-i spin liquid #19: up to 2nd n.n. bonds needed 

The mean-field ansatz {uij} for Z2 spin liquid t^17 is 
summarized in TABLE [III up to 2nd n.n. bonds. The 
corresponding spin liquid has a Z2 gauge structure, if 
and only if up = (3it^ +1^2^'^ with /?i, /32 7^ 0, so that the 
IGG of this mean-field ansatz is a Z2 group {±r°}. This 
is the only gapped Z2 spin liquid near the u-RVB state, 
that can be realized in a mean-field ansatz up to 3rd n.n. 
bonds. 

It's straightforward to check that 2nd n.n. bond up — 
PiT^ + /32T^ open up a mass gap M ^ /32/^'^ (E) u^ (E) t^ ~ 

P2Mf. 



corresponding spin liquid has a Z2 gauge structure, if and 
only if [up,ug] ^ 0, so that the IGG of this mean-field 
ansatz is a Z2 group {±r°}. wg = (1, 2, 0) is the 9th n.n. 
mean-field bond. In this Z2 spin liquid, the symmetry- 
allowed consistent mean-field bonds for 6th, 7th and 8th 



n.n. are: 



UgSE (2,0,0) (xt2, 
UT = (2,0,1) ex ir°, 
M8=(0,2,l)(x irO. 



d. Z2 spin liquid #22: up to 9th n.n. bonds needed 



The mean- field ansatz {wij} for Z2 spin liquid #17 
is summarized in TABLE |lTl up to 9th n.n. bonds. The 



It's straightforward to check that 2nd n.n. bond up 



ISt"^ open up a mass gap M ~ /3/x^ (g) z/^ (g) r^ 



PMf. 
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